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In this paper, we derive a Darboux transformation of the Hirota and the Maxwell-Bloch (H- 
MB) system which is governed by femtosecond pulse propagation through an erbium doped fibre 
and further generalize it to the matrix form of the n-fold Darboux transformation of this system. 
This n-fold Darboux transformation implies the determinant representation of n-th new solutions 
of {E^^\p^^\ri^^^) generated from known solution of {E,p,ri). The determinant representation of 
{E^^\p^^\ri^^^) provides soliton solutions, positon solutions, and breather solutions (both bright 
and dark breathers) of the H-MB system. From the breather solutions, we also construct bright 
and dark rogue wave solutions for the H-MB system, which is currently one of the hottest topics in 
mathematics and physics. Surprisingly, the rogue wave solution for pandrj has two peaks because 
of the order of the numerator and denominator of them. Meanwhile, after fixing time and spatial 
parameters and changing other two unknown parameters a and /3, we generate a rogue wave shape 
for the first time. 
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I. INTRODUCTION 

In the past four decades, nonlinear science has ex- 
perienced an explosive growth with the invention of 
several exciting and fascinating new concepts such as 
solitons, dromions, positons, rogue waves, similari- 
tons, supercontinuum generation, complete integrabil- 
ity, fractals, chaos etc. Many of the completely inte- 
grable nonlinear partial differential systems (NPDEs) 
admit one of the most striking aspects of nonlinear 
phenomena, described as the soliton, a universal char- 
acter and is of great mathematical interest. The study 
of the solitons and other related solutions like positons 
have become one of the most exciting and extremely 
active areas of research in the field of nonlinear sci- 
ences. 

Among all concepts, in addition to solitons and 
positons [1-4 , rogue waves have also been not only 
the subject of intensive research in oceanography [5- 
|7] but also they have been studied extensively in sev- 
eral other areas, such as matter rogue wave ^ in 
Bose-Einstein condensates, rogue waves in surface and 
space plasmas [10], financial rogue waves describing 
the possible physical mechanisms in financial markets 
and related fields [11]. In some of the above fields, 
soliton system such as nonlinear Schrodinger (NLS) 
equation [12 , derivative NLS system [13, 14 and so 
on are considered and reported to admit rogue wave 
solutions under a certain specific choice of parame- 



ters. It has been proved that modulational instability 
is one of the main generating mechanisms for the rogue 
waves [12 -17 and can be well-described by the analyt- 
ical expressions for the spectra of breather solutions 
at the point of extreme compression. 

In 1967, McCall and Hahn [18 explored a special 
type of lossless pulse propagation in two- level resonant 
media. They have discovered the self-induced trans- 
parency (SIT) effect which can be explained by using 
the Maxwell-Bloch (MB) system. If we consider these 
effects in erbium doped nonlinear fibre, the system 
will be governed by the coupled system of the NLS 
and the MB equation (NLS-MB svstem) [19^23) . 

Rogue waves have been reported in different 
branches of physics, where the system dynamics is 
governed mostly by a single nonlinear partial differ- 
ential equation [l3l [TJ . But our main interest is to 
analyze the possibility of rogue waves in coupled non- 
linear systems. The higher-order NLS and Maxwell- 
Bloch (HNLS-MB) system as a higher-order correc- 
tion of NLS-MB system were shown to admit Lax pair 
and soliton- type pulse propagation [25H27] . Kodama 
[24] has shown that with a suitable transformation, 
higher order NLS equation can be reduced to the Hi- 
rota equation [28 whose rogue wave solution has al- 
ready been reported in [29l [S^ . In a similar way, after 
suitable choice of self-steepening and self-frequency ef- 
fects, we obtain the H-MB system in the following 
form [31 : 

= ia{^Ett + \E\^E) + l3{Ettt + Q\E\^Et)+2p, 

(1.1) 

Pt = 2iujp + 2Eri, (1.2) 
r]t = -{Ep* + E*p), (1.3) 
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where E is the normalized slowly varying amplitude 
of the complex field envelope, p is the polarization. 



T] means the population inversion, (co',a,/3) are three 
real constants and * represents complex conjugate. /3 
represents the strength of the higher order linear and 
nonlinear effects. 

The H-MB system has been shown to be integrable 
and also admits a Lax pair and other required prop- 
erties of complete integrability [31 . Among many 
analytical methods, it is well known that the Dar- 
boux transformation is one of the efficient methods to 
generate the soliton solutions for integrable systems 
[32]. The determinant representation of n-fold Dar- 
boux transformation of the Ablowitz-Kaup-Newell- 
Segur (AKNS) system was given in [33l[34]. The main 
task of this paper will be to construct n-fold Darboux 
transformation of the H-MB system and find differ- 
ent kinds of solutions of the H-MB system using the 
Darboux transformation. 

The paper is organized as follows. In section 2, the 
Lax representation of H-MB system is introduced. In 
section 3, we derived the one- fold Darboux transfor- 
mation of the H-MB system. In section 4, the gen- 
eralization of one-fold Darboux transformation to n- 
fold Darboux transformation of the H-MB system will 



be given. Using these Darboux transformations, one 
soliton, two soliton and positon solutions are derived 
in section 5 & 6 by assuming trivial seed solutions. 
In section 7, starting from a periodic seed solution, 
breather solution of the H-MB system is provided. A 
Taylor expansion from breather solution will help us 
to construct the rogue wave solution in section 8. Sec- 
tion 9 is devoted to conclusion and discussions. 



II. LAX REPRESENTATION OF THE H-MB 
SYSTEM 

In this section, we will concentrate on the linear 
eigenvalue problem of the Hirota and the Maxwell- 
Bloch(H-MB) system. The linear eigenvalue problem 
is expressed in the form of the Lax pair U and V as 



= U^, = (2.1) 
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is an eigenfunction associated with eigenvalue param- 
eter A of the linear Eq. (2.1 ), and Vi denotes the coeffi- 
cient matrix of term A*. We obtain the classical Hirota 
and the Maxwell-Bloch system when a = 2,/3 = —1. 
Being different from AKNS system, only a part of the 
V matrix is polynomials in terms of E and its t deriva- 
tives in this system. Using the above linear system of 
the H-MB system, one-fold Darboux transformation 
will be introduced in the next section. 



HI. ONE-FOLD DARBOUX 
TRANSFORMATION FOR THE H-MB 
SYSTEM 



In this section, we construct and prove the one-fold 
Darboux transformation for the H-MB system. First, 
we consider the transformation about linear function 
$ in the form 



(3.1) 



where 
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The new function satisfies 



A2: 



(3.3) 
(3.4) 



Then the matrix T should satisfy the following iden- 
tities 



Tt + TU 

n + TV 



U'T, 
V'T. 



(3.5) 
(3.6) 



Substituting the matrices A and S into Eq. (3.5) and 
comparing the coefficients of both sides will lead to 
the following conditions 



^12 CL21 0, (aii)t = {a22)t 0. 



(3.7) 



For our further discussions, we choose A = I and 
T = {XI — S). The relation between old solutions 
{E^p^T]) and new solutions {E'^p'^r]')^ which is called 
Darboux transformation, can be obtained by using 
Eqs. (|3^ an d K^ . 

FroiXL Eq. (3^), we have 



E' E-2is 



12, S21 — —Si2 



(3.8) 



St = 



E 

-E* 



s-s 



(3.9) 



Similarly, using Eq. (3.6), we obtain the following set 
of relations 



-S, + {XI - S){X^Vs + AV2 + XVi + Vo 



= (AV3 + + w; + 



A + w 



){\I-S). 



(3.10) 



Multiplying both sides of Eq.( |3.10| ) by A/ - 5 will lead 
to 



-S,{\ + oj) + (AJ - 5)(-4iA3(A + w)(T3 
+A2(A + uj)V2 + A(A + uj)Vi + l^o(A + w) + iV-i) 
= {-Ai\^{\ + uj)az + \^{\ + oj)V2 

+A(A + uj)Vi + Fo'(A + w) + iVLi){\I - S). 

Collecting the different powers of A, we obtain the 
following set of identities 
AO: 



S, = {V^+iuj-WLi)S-S{iuj-^V-i + Vo), 



(3.11) 



A: 



= {ojVo + iV-i)-S{ojVi+Vo) + {LoVi + V^)S 
+{-LoV^-iVLi). (3.12) 



A3: 



{ujVi +Vo)- S{ujV2 + Vi) (3.13) 
[ioVi + T^o') - i^V^ + Vi)S. 



{loV2+Vi)-S{loVz + V2) (3.14) 
{u:Vi + Vi) - {uV^ + Vi)S, 



Vi = V2-[S,V3\. 



(3.15) 



From the above identities, after simplifications, we 



get 



E' = E-2is 



12, 



VU = {S ^ u)V-i{S ^ 



(3.16) 



(3.17) 



which gives one-fold Darboux transformation of the 
H-MB system later. 
We suppose 



S = HAH-^ 



(3.18) 



where A 



H 



Ai 
A2^ 

$i(Ai,t,z) ^i(A2,t,z)\ ^1,2' 

<l>2(Ai,t,z) <l>2(A2,t,z)y • \^<I>2,1 ^2,2^ 

In order to satisfy the constraints of S and VLi 
which is similar to V-i, i.e. S21 = —8^2^ following 
constraints will be used 



A2 — A^, Sii — 5^2, 

/$i(Ai,t,^) -$5(Ai,t,z; 



H 



^2(Ai,t,z) ^l{Xi,t,z) 



(3.19) 
(3.20) 



Aft er t edious calculations, Eqs. ( |3.16| )-( |3^ ) and 
Eq. ([2I]) win lead to E q. (|3^ and Eq. 13^ , i.e. 
the transformation Eq. ( |3.16 ) and Eq. (3.17) with 
the conditions Eq. (3.19), Eq. (3.20) is the Darboux 



transformation of the H-MB system. 

The detailed form of one-fold Darboux transforma- 
tion of the H-MB system in terms of eigenfunctions 
will be given in the next section. 



IV. DETERMINANT REPRESENTATION OF 
n-FOLD DARBOUX TRANSFORMATION 

In this section, we will construct the determinant 
representation of the n-fold Darboux transformation 
of the H-MB system. For this purpose, we introduce 
n eigenfunctions 



*2,i 



A=Ai> 



1,2,. ..,2n 



(4.1) 



3 



with the following constraint on the eigenvalues 
A2n-i = ^2n ^^'^ ^he reduction conditions on eigen- 

functions as $2,2n = *l,2n-l«"'^ $2,2n-l = -*l,2n- 

For our further discussions, this reduction condition 
has been used. 

For completeness, as the simplest Darboux trans- 
formation, the determinant representation of one-fold 
Darboux transformation of the H-MB system will be 
introduced in th e following theorem using identities 
dale] ) and ( |3l7| ). 

The one-fold Darboux transformation of the H-MB 
system is expressed as 



where 
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(At - Ai)$i,i$;,i - p*{X* - Ai)2$2 

+p[{-OJ - Ai)|$i,i|2 + {-OJ - At)|$2,llT), 

(4.3) 



.[1] = 



[r,([(c^ + Ai)|$i,i|2 + (w + At)|$2,il'] 



[(c^ + At)|*i,i|2 + (a; + Ai)|$2,il'] 
+(AI-Ai)2|$i,i|2|$2^|2) 

-p*(Ai-At)$i,i$;,i[(w + At)|$i,i|2 

+(w + Ai)|$2,iP]+H(w + Ai)|*i,iP 

+(w + ADI*2,il'](Ai-At)$ti$2,i], 



Ta, = [(w + Ai)|$i,i|2 + (t^ + AI)|$2,l|'] 
[(W + At)|$i,i|2 + (C. + Ai)|$2,ll'] 

-(At-Ai)2|$i,i|2|$2^i|2. (4.5) 



It can be easily proved that the new solution rj^^'^ 
is always real. This one-fold transformation will be 
used to generate the one-soliton solution from trivial 
seed solutions of the H-MB system. Also this one-fold 
Darboux transformation can be further generalized to 
construct the n-fold Darboux transformation of the 
H-MB system which is proposed in the following the- 



Ai)|$i,ir + (-w-AI)|*2,il ] orem. 



Theorem IV. 1. The n-fold Darboux transformation 
of the H-MB system can be represented as 



?n(A; Ai, A2, A3, A4, . . . , A2„) (4.6) 

= A"/ + 4"1iA"-i + --- + 4"'a + 4"' (4.7) 
^ J_ AT„)ii (T„)i2 
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The following identity can be found 



Tn(A; Ai, A2, A3, A4, . . . , A2n)|A=Ai 



0, (4.9) 



where i = 1,2, 3,..., 2n. Similarly, the Darboux trans- 
formation for T]) can be constructed by using the 
following identities 



T„t+T„U = U^^'^T^, 
which can be further simplified as 



tn|A= 



X=—oj ' 



(4.10) 
(4.11) 



(4.12) 



(4.13) 



The proof for Eq. ( |4.13 ) is quite complicated in 
general when compared to the proof of one and two 
fold Darboux transformation. However, this will be- 
come simple if we know the origin of [33]. Be- 
cause if we treat n-fold Darboux transformation as a 
generalization of (n — l)-Darboux transformation, the 
transformation will be the multiplications of n one- 
fold Darboux transformations at A = —uj. It is easy 



to prove that these multiplications can have a deter- 
minant representation as mentioned above. 

The n-th new solution after the n-fold Darboux 
transformation of the H-MB system will be 



£["1 = i? + 2i(41i)i2, (4.14) 

pN = (2r?(T„)ii(T„)i2-p*(T„)i2(T„)i2 

+p(T„)n(T„)n) (4.15) 

/((r„)ii(T„)22 - (T„)i2(T„)2l)| A=— a; 7 

7/["i = (r?((T„)n(r„)22 + (r„)i2(r„)2i) (4.16) 

-p*(T„)i2(T„)22 +p(T„)n(r„)2i) 

/((T„)ii(T„)22 - (T„)i2(T„)2l)| A=— a; 7 



where (t\^}_i)-^^ is the element at the first row and 

second column in the matrix tj^l^. So far, we dis- 
cussed about the determinant construction of n-th 
Darboux transformation of the H-MB system. As an 
application of these transformations, soliton and posi- 
ton solutions of the H-MB system will be constructed 
in the next section. 
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V. SOLITON SOLUTIONS OF THE THE 
H-MB SYSTEM 

In this section, having obtained the Daurboux 
transformation for our system, our next aim is to con- 
struct the one sohton solution of the H-MB system by 
assuming suitable seed solutions. We assume trivial 
seed solutions as E = 0,p = 0,7^ = 1, then the linear 
system becomes 



where 



$ = 
U = 



$1 

$2 

-iX 
iX 



'ii/SX^ - aiX^ 
^ ~ ' -4/?iA3+mA2 



i \ 



(5.1) 
(5.2) 



(5.3) 
(5.4) 

(5.5) 



From the above system, we construct the explicit 
eigenfunctions in the form 



= g-iAt+(4^iA3-aiA2 + ^)2;+ 
^ = g*At+(-4/3iA2+aiA2-^)2- 



_ XQ+iyQ 



where xq, Vq and are all arbitrarily fixed real con- 
stants. Substituting these two eigenfunctions into the 
one-fold Darboux transformation Eg. (4.2), Eq. (4.3), 



Eq. (4.4) and choosing A = ai + xq = 0, ?/o = 



0, ^ = 0, then the following form of soliton solutions 
are obtained: 



E = 2pie-^'Bsech{2pi-), 



(5.6) 



iPiiiai - + w)e-2s + (aj + l3ii + uj)e-'^B) 



B 



sech\2l3i-), 



V = l-^sech\2/3^^), 



(5.7) 
(5.8) 



where A^B^C^D^F are explicitly given in Appendix 
I. 

If we choose a = 1, /3 = 0, the one soliton solution 
is just the soliton solution of Eq. (15) of NLS-MB sys- 
tem mentioned in [22^ with ai = — pi, Pi = vi. Sim- 
ilarly, substituting these two eigenfunctions into the 



one -fold Darboux transformation eq.(4.2), eq.(4.3), 
eq.(4.4) and choosing a = 2, /3 = —1, then the one- 



solition solutions of the classical H-MB system can 
be obtained whose evolution is shown in Figjl] which 
clearly indicates that E and p are bright solitons be- 
cause their waves are under the flat non-vanishing 
plane whereas is a dark soliton. 

Now let us discuss the construction of the two- 
soliton solutions of the H-MB system. For this pur- 
pose, we have to use two spectral parameters Ai = 
ai + iPi and A2 = 0^2 + ^^2- After the second Dar- 
boux transformation, we can construct the two so- 
lition solutions. As the general form of two soliton 
solution is quite tedious in nature, for simplicity, we 
are giving only the two soliton solution of E with 
u = 1.5, Qfi = 0.5, /3i = 1,0^2 = l,/^2 = 1.5, a = 2 
and /3 = — 1, 



E2-S0I = 

_^j^-|^yg-0.01176470588i(85t-12582+1815i2+255it) 
_(^35 ^ -|^g^jg-0.01176470588i(170t-43852+204i2+170i^) 
_Yg0.01176470588i(-85t+12582+1815i2+255it) 

-\-{18i — 12)6^-^-^-^'^^^'^^^^^*'^"-^'^^^+^^^^^+^^^*^+-^'^^*^^ 

_2g^g0.01176470588i(-85t+12582+1815i2+255i^) 
+20ie~^"^"^"^'^^^'^^^^^*'^^^^~"^^^^^^"^^"^^*^^^^^*^^) 

/(-24co5(0.9999999998t - 36.78823529z) 
^26cosh{t + 18.95294118^) 
^2cosh{-5t - 23.75294118z). 

We also constructed the two soliton solution for p 
and 77 in a similar manner. For completeness, instead 
of giving complicated forms of p and 77, the graphical 
representation of them is shown in Figj2] 



VI. BRIGHT AND DARK POSITON 
SOLUTIONS OF THE H-MB SYSTEM 



In the case of two soliton solution constructed 
above, if the second spectral parameter A2 is assumed 
to be close to the first spectral parameter Ai, and do- 
ing the Taylor expansion of wave function to first order 
up to Ai will lead to a new kind of solution which is 
called as degenerate soliton [34 - smooth positon so- 
lution. "Positon" was coined by Matvee [1, 2, 4 for 
the Korteweg-de Vries (KdV) equation by the same 
limiting approach. Note that the positon of the KdV 
is a singular solution. For the construction of positon 
solutions, following four linear functions are needed to 
construct the second Darboux transformation and to 
generate the positon solutions. 
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^1,1 

^2,1 
^1,3 
^2,3 



= e 



iAit+(4/3iA?-mA2 + ^^)2H 



= e 



iAit+(-4/3iA?+aiA?-^^)2-^a±^+i6' 
^-iA3t+(4/3iA^-mA^ + ^)^+^o+lm ^ 

A3t+(-4/3iA^+mA2- ^o±^+i^ 



i(Plz + Qlt) 



2d' 



(7.4) 

R)eMb-2u.)(x+u.) ^ (7.5) 



Now we take A3 = Ai + e(l + i) and use the Tay- 
lor expansion of wave function ^3 and ^4 up to first 
order of e in terms of Ai. For example, choosing 
uj = 1.5, ai = 0.5, /3i = l,a = 2 and j3 = —1, the 
positon solutions Ep is constructed in the form 



{15S2zcosh{2t + 2A0z) - 6SSizsinh{2t + 2Az) 
-50icosh{2t + 2.4^) + 100itsinh{2t + 2.4^))/ 
(400t^ + 119296^^ - 5504tz + 50 + 50cosh{4.Sz + 4t)). 

In this case, the pictorial representation of the posi- 
ton solutions {Ep^pp^ rjp) of the H-MB system is shown 
in FigjS] 

From the above figures, we observe that two peaks 
of the positon solutions are at same height which is 
different from the two-soliton solution. Meanwhile the 
two waves depart at a relatively less speed after their 
collision. This is also different from the two-solitons, 
which depart at a fixed speed. Here again, we find 
that E and p are bright positon solutions whereas r] 
is a dark positon in all three cases discussed above. 



VII. BRIGHT AND DARK BREATHER 
SOLUTIONS OF THE H-MB SYSTEM 

In the last two sections, soliton solutions and posi- 
ton solutions have been generated for the H-MB sys- 
tem. In this section, we will now focus on a new and 
different kind of solution which is also derived from 
periodic solutions through Darboux transformation. 
The resulting periodic solutions can be called breather 
solutions. Now, let us assume the seed solutions as 
E — de^P^p = ifE^T] = 1^ p = az -\- bt^ which admits 
the constraint in the form 



where Pi,Qi,P2,Q2 are polynomials independent of 
t and z and their complete expressions are given in 
Appendix II. 

If we use the above two wave functions to con- 
struct the two new functions hi and /12 as in [22. , then 
the resulting new solutions obtained through Darboux 
transformation are found to have no meaning. There- 
fore, we would like to construct more complicated but 
physically meaningful solutions in the following part. 
By combining these two wave functions, we derive the 
new functions hi and /12 as follows 

hi := 51(A) -<7i(A*)*,/i2 :=<72(A) + <?i(A*)*, (7.6) 

where 

gi := MR) + M-R).92 := MR) + M'R)- (7.7) 

It can be proved that hi and /12 are also the solu- 
tions of the Lax equation with A := ai + i^i. Using 
these two wave functions hi and /12 in the one-fold 
Darboux transformation will lead to the construction 
of breather solutions of the H-MB system. To sim- 
plify the calculations, we take b = —2ai and use the 
second Darboux transformation discussed in the last 
section, then the final form of the breather solution 
E}) is obtained in the form 



Eb = (d'(e^^ +e^^^) + 2d/3i(e^^ +e^^) 



-dMe^' + e^') + 2/3i(^ - f3i)e^' 

-2/3i(^ + /3i)e^0/ 

G H 
{2dcosh{2iwPiZY) - 2Picosh(2w^)), 



(7i 



where 



(7.9) 



2a + ab"^ - 2a(f + 2pb^ 



Defining 



- 12P(fb - 4:f = 0, (7.1) 
/6-2w/ + 2 = 0. (7.2) 



R := A/-52 -45A-4A2 -4^2, 



(7.3) 



the following wave function <I> 
terms of R as 



/ I is obtained in 

^2; 



X = (ai + Pii + uj){-ai + Pii - uj), (7.10) 

r = (ai + uj){ai + /3ii + uj){-ai + /3ii - cj),(7.11) 

and A}){w),Cb^ Dij, F}){w), and Gb^Hij are polynomi- 
als of t, 2:, cj, a, 6, d, a, /3, ai, Pi which are defined in 
Appendix-II. 

Similarly, the breather form of p and rj can be con- 
structed. For example, after taking values 

UJ = 0.5, d = 0.5, ai = = l,a = 2,(3 = -1, 

(7.12) 
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the breather solution of the H-MB system is plotted 
m Fig(4]and FigjS] 

Similarly, for the next choice of parameters uo = 
1.5, ai — 0.5, /3i = 1, a = 0, /3 = 1, the breather 
solutions of the complex modified Korteweg-de Vries 
(CMKdV)-MB system are obtained. The picture of 
breather solutions is shown in Fig|6] In addition to 
the above observation, we also find how the values 
of (a, P) changes the direction of breather solution 
T] in the (t, z) plane, see Figj?} To the best of our 
knowledge, we observe this effect for the first time. 

Having constructed bright breathers for E and p 
and a dark breather for 77, in the next section, our aim 
is to discuss the construction of rogue wave solutions 
of the H-MB system which is in fact one single period 
of breather solutions. 



Pr 



8(1 - 3iz) 



0.5i(-52+4t) 



, (8.2) 



M + (2t + 17z)2 + 36^2 2 
g0.5i(-5^+4t)^^Q3^4g3Qy3-^4 ^ io3i2500zz^ 

-5312500^ - 625000t - 152343760^^ 

-1875000t^z - 31875000tz^ 

+70546875zt^z^ + 5312500it^z + 431640580itz^ 

+1875000itz + 156250^^"^ - 1562500/(156250^"^ 

+25390625^^ + 312500t^ + 431640580tz^ 

+70546900t^z^ + 5312500t^z 

+156250 + 5312500tz + 1031493636^"^), (8.3) 

(16t^ + 64^2 + 722At^z^ + 105625^"^ + 2992^^ 

+44200tz^ - 16 + 544t^z + 896zt) 

7(16^"^ + + 7224t^^^ + 105625^"^ 

+2600^^ + 44200t^^ + 16 + bUt^z + 544zt). 

(8.4) 



VIII. BRIGHT AND DARK ROGUE WAVES 
IN THE H-MB SYSTEM 



In this section, using the limit method of the NLS 
equation, we construct the rogue wave solutions of the 
H-MB system [34] . This kind of solution only appears 
in some special regions of time and distance and then 
will be drowned in one fixed non- vanishing plane. If 
we d o the Taylor expansion to the breather solution 
(7.8) around j3i = one rogue wave solution of E 



will be obtained and rogue waves for p and rj can also 
be constructed in a similar way. In the following, for 
brevity, we only report the rogue wave for Er in the 
form 



Er 



exp(- 



ai 



- uo 



-(Szpa\^Szpaluj 



A{t,z)d 
B{t,z) 

—2zaa\uj + 2z — 2aituj + zadPuj — 12zPdPaiUJ 
-2zaal - 2alt + zad^ai - I2zpd?al)), (8.1) 



From Eq. ( 8.2 ) , it is clearly observed that the height 
of the background of l^^p is and the orders of the nu- 
merators and denominators of Pr and rjr are four. Be- 
cause of these reasons, the graphs of the rogue waves 
for Pr and rjr have double peaks which are shown in 
FigjS] and the corresponding density graph is plotted 
in Fig(9| 

For further understanding of our observations, we 
enlarge the above density in Figj9j some zoomed por- 



tions of the above figures are clearly shown in Fig 10 
From the graph of \p\^ shown in Fig 10, we find one 



cave appears on top of the single peak with two caves 
on both sides of the peak. 

To realize the significance of the different parame- 
ters a and /3, we also consider the case, when a = 

0, /3 = 1, i.e., the CMKdV-MB system, see Fig{Tl] 
From the figures, we also find that the parameters a 
and 13 will change the shape, pulse width, etc., of the 
rogue wave. Therefore, in the following, we fix time t 
and distance z to see the role of parameters a and P 
and their impact on rogue wave dynamics. 

We keep a and P as arbitrary parameters and 
choose uj = 0.5, d = 0.5, ai = —1, /3i = 1, t = 1, z = 

1. This is also a rogue wave whose graph is portrayed 
in Fig 12 Here we provide only the specific form of 
the rogue wave solution Erai3 as 



where A{t^ z)^ B{t^ z) are polynomials of 
t, a, 6, d, a, /3, Qfi, /3i which are defined in 

Appendix-IIL 

When we take u = 0.5, d = 0.5, ai = —1, /3i = 
l,a = 2,/3 = —1, the final form of the rogue wave 
solutions will be 



Era^ = e-^^^+20^+^«)(40a- 17^2 + 156/3 -40 
-585/3^ - 192a/3 + i{32 + 8a + 96/3))/ 
(1170/3^ - SOa + 34a^ - 312/3 + 112 + 384a/3). 

(8.5) 

This implies that after fixing the values of t and 
the solution depending on parameters a, /3 is also in 
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the form of a rogue wave, which is observed for the first 
time. This will give us some idea about how to modify 
the parameters a and {3 to visualize our theoretical 
results in terms of the experimental results in optics. 
From the above, one can easily conclude that E and 
p are bright rogue waves and is a dark rogue wave. 

All the solutions mentioned above including posi- 
tons and rogue wave solutions are indeed solutions of 
the H-MB system which are verified by using the sym- 
bolic computation software MAPLE. 



IX. CONCLUSION AND DISCUSSIONS 

In this paper, after a suitable choice of self- 
steepening and self-frequency shift effects, we have 
derived the Darboux transformation of the H-MB 
system which is governed by ultra-short pulse prop- 
agation through an erbium doped nonlinear optical 
waveguide and further generalized it to the matrix 
form of an n-fold Darboux transformation, which im- 
plies the determinant representation of [E^'^'^ , p^'^'^ , /^["^^ ) 
generated from the known solution {E^p^rf). By 
choosing some special eigenvalues \2n-1 = ^2n 
and eigenfunctions using the reduction conditions 



^l,2n-l = ^2,2n. ^2,2n-l 



"^1 2n5 determinant 



representation of [E^'^\p^'^\r]^'^^) provided some new 
solutions of the H-MB system. As examples, soliton 
solutions, breather solutions, and rogue wave solutions 
of the H-MB system have been constructed explicitly 
by using the Darboux transformation from trivial and 
periodic seed solutions. The rogue waves show inter- 



esting characteristics which might attract physicists 
to observe them in experiments with higher order op- 
tical effects in the femtosecond regime. The interest- 
ing characteristics obtained contain the following two 
sides: (i) The rogue wave solution for p, 77 is surpris- 
ingly found by us to have two peaks because the order 
of the numerator and denominator of p, r] in eq.(8.3) 



and eq.(8.4) is four and (ii) after fixing the time and 



spatial parameter and by changing other two unknown 
parameters a and /3, we fi nd a rogue wave shape also 
arises out as shown in eq.(8.5). This is the first time 



that this phenomenon is obtained to the best of our 
knowledge. Still, there are a few interesting questions 
which are still unclear. For example, the physical in- 
terpretations and observation of higher-order positon 
solutions, the role of higher-order rogue waves solu- 
tions and their applications in physics, in particular, 
the connection between rogue wave solutions and su- 
percontinuum generation through modulation insta- 
bility or soliton fission, etc., with higher-order optical 
effects. 
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X. APPENDICES 



Appendix I: In this appendix, we are providing explicit expression for A, B, C, D, and F 



-UzPaluj'^ + z + SzPaiPfuj + taj + tp^ + tuj'^ + 2taiuj - I2zpa\ + 2zaal + 4z/3/3^, 
B = al^2aiuj ^uj'^ ^ Pl, 

C = 2ta\uj + taiuj'^ + t^lai - Azpa\ + zaa\ - za^f - zai + taf 

^24z(3al(3luj + UzpaifSfuj'^ - 2zaai(3luj - 8z(3ajuj - Az^Safuj'^ 
-\-Szf3alPl + 12z/3P^ai + 2zaaluj + zaajuj'^ - zapfuj'^ - zoo, 



D := -24PizPaluj + APizaajuj + 2PizaaiUJ^ - 12 (3 iz (3 ajuj'^ + SPfz(3aiUJ + (3iz + (3ft + 4pfz(3 
+/3ituj'^ + /3ital - 8f3fzpal + 2f3fzaai + 4f3fzpuj'^ + 2pitaiuj - 12f3izpat + 2f3izaal 
-\-i{tplai + zaaf — zol\ — zuj ^ ta\ + taiuj'^ + 2ta\ijj — zaPfuj'^ 
-4z(3aluj'^ - 8z(3ajuj - 4z(3al - za(3f - 2zaai(3luj 

-^SzPalPl + UzPPfai + 2zaaluj + 24zpalpluj + I2zpaipluj'^ + zaaluj'^), 
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= 24(3izl3aluj - A/Sizaaluj - 2l3izaaiuj^ + UPiz/Saluj^ - SPlz/SaiUJ - Piz - /Sft - 4/3^ 
-Pitoj'^ - Pitaj + SPfzPaj - 2Plzaai - Apfz^uj'^ - 2Pitai(jj + UPiz^aj - 2Pizaa' 
-\-i{tl3fai + zaa^ — zol\ — zuj -\- ta\ + taiuj'^ + 2ta\uj — za^fuj'^ 

—4:z(3af(jj^ — Sz(3afuj — Azfia^ — za(3f — 2zaaiP\uj 

-^SzPalPl + UzPPfai + 2zaaluj + 2AzpalPluj + I2zl3aipluj^ + zaaluj'^), 
Appendix II: In this appendix, we are furnishing the expression for Pi, Qi, P2, Q2 

Pi = l2l5dH^\-Aad?\uj + AI5h^\uj + 2ad?huj + 2ad?h\ 

^2ah^\uj - 24P(fb(jj'^ + UpdH'^oo -SX-Soo- Aad^uj'^ - 2/36^A 
-2l3y^uj + 2ab^uj'^ + 4^&^a;^ - ab^X - ab^oj - 24pdHX(jj + 2iujPRb^ 
-^2iXpRb^ + XaRb^i + ujaRb^i + WiX'^uj^Rb + 4iR + Si^Xbuj'^R 
-SipXb'^ujR - Aid^puoRb - Aid^pXRb - AiabXRuj + Sid'^pXRuj 
-Aipb'^uo'^R - 16i/3X^Ruj - 2iaX^Rb - 2iabuj'^R - mpX^co'^R 
-AipxH'^R + UaX^Ruj + SipX^Rb + m^Puj'^R + AiaXoo'^R, 

Qi = 2b^uj - 4bX(jj - 4b(jj'^ + 2b'^X - 2iRbuj + URujX + URuj'^ - 2iRbX, 

P2 = -12pd'^b'^X^4ad'^Xuj-4l3b^Xuj-2ad'^buj-2ad'^bX-2ab'^Xuj^24f3d'^buj'^ 
-I20b^uj + 8A + 8cj + 4.ad^uj'^ + 2^6^A + 2pb^uj - 2ab^J^ - 4.pb^uP 
-\-ab^X + ab^uj + 24l3d'^bXuj + 2iujl3Rb^ + 2iXl3Rb^ + Acei^fe^i + (jcei^fe^i 
+16zA^cj/3i?6 + 4ii? + 8i(3Xbu'^R - 8i(3Xb'^ujR - Aid^pujRb - Aid^pXRb 
-AiabXRcj + 8id^f3XRuj - 4if3b^uj^R - 16if3X^Ruj - 2iaX^Rb - 2iabuj^R 
-WipX^uj^R - Ai^X^b^R + AiaX^Ruj + Si/SX^Rb + 8id^l3uj^R + AiaXuj^R, 

Q2 = 4.bXuj + Abuj'^ - 2b^(jj - 2b^X - 2iRbuj + AiRcuX + AiRcu'^ - 2iRbX, 

Appendix III: In this appendix, we are providing the expression for 
Ab{w), Cb, Db, Fb{w), Gb Hb,A{t,z)andB{t,z) 

= —2wl3iz — 12zPa\ujwPi — 72 z (3 aluj'^w Pi + ^zaa\ujw Pi + ^zaaiUj'^wPi 

-24:zPujwaiPf - 2AzpaiUJ^wPi + 2zaa\wPi - 2AzPa\wPi + 2zauo^wpi - 24:Zpwa\ 
-\-2zawaiPf + 2zawujpf - {azaf + azuj^ + bta^ + btuj^ + Sazajuj 
-\-3azaiu;'^ + Sbtajoj + S&tceio;^ + azaiPf + a^^cj/^^ + btaiPf + btouPl), 

Cb = -2zuj'^ - 2zpl + 2a^t - 2za\ - Azujai + Qaltuj + 2^aaf - 8^/3a^ 
-\-Qaituj + 2ceitct; — zad uo — zad + QzaaiCJ + G^ao^io; 
+2^aa?cj^ - 24zPaluj - 2AzPa\uj'^ - Sz^afuj^ + I2zpd^a\ 
—Szad^uj^ai — 3zad^aiUJ + SGz^d'^a^uj + 36z/3(i^a^cj^ 
^UzPd'^aiu^ - Szpajpf + 2zaa?/3f + 2aitujpl + 2a?t/3f 
^Uzd^pwujpl + Aizaaiwujpl + Aizd^ pwaipl + 24i2;/3aicj^it;/3i 
+12izd^/3cc;^'L<;cei + I2izd^ Pwa\uj + 2itwujPl + 2itwaiPl 
-\-^itw(jj^ai + ^itwa\uj + Aizaa\w — 2AizPa\w + 12zPd^a\ 
P1 + 2zaa^a;/3^ — zad^ai/^^ — zad^ujp\ — Sz^afuj^i 
^I2zpd^aiujpl - 12izpa\ujw - 12izpalu'^w - 2Uzpalu^w 
-\-4:izawalPf + 12i2;aafcji(; + 12izaa^co'^w; + 4izaaicj^i^ 
+4izc^^^cc;^'L<; + Uzd^ Pwal - IQizpalwpl + Sizpuj^wPl + Siz^wai 
Pf + SizPwouPi + 2i2:i(;cei + 2izw(jj + 2itwuj^ + 2itwai, 
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-\-Siz/3cu'^w/3f + 2i2:i(; + Aizaa\w + Aitwaiuo — 4:Siz/3aiCuw — 2AizPa\(jj^w 
-\-2itw(jU + 2itwai — 2MzPa-^w + Si^^c/ /Scuwai — aza^ — azou — bta^ 
—btu^ — 2azaiuj — 2htaiuj — az^f — ht^l + 8izPw/3f + 16iz/3ujwai/3f 
^Aizawaipl - l^izpalwpl + Uzd^(3w(3f + 2it^(;/3^, 
Fb = -2zuj'^ - 2zPl - 2wPiz + 2aft - 2zot\ - Azuai + Galtuj + 2za 

— 8z/3ai + ^a\tuP + 2aitco^ — zad^co^ — zad^a\ + 6zaafco + 62:0; 
o;^ + 2zceQf^a;^ - 24zPaluj - 24zf3ajuj'^ - Sz^a^uj^ + Uz^d'^ 
ol\ - ?>zoLd^iJ'OLx - Szad^ajuj + 36zpd^aluj + 36z(3d^aluj^ + 12z/3c^^ai 
cj^ — 72zPalujwPi — 72z(3aluj'^w(3i + 6zaalujw(3i + Gzaaicj^i^/^i 
-24:zPuwaiPf - 24:zPaiUJ^wPi + 2zaa\wPi - 2AzPa\wpi + 22;acj^ 

- 24:zPwalPf + 22:aw;ai/3i + 22;ai(; 
cj^^f - SzPafPl + 2zaalPl + 2aita;/3^ + 2a^t/3^ + 12z/3(i^a^/3^ + 2z 
aajuoPf - zad^aiPl - zad^uj^l - Sz^aluj^f + Uz^d'^aiuj^f, 



Gh = -l-12Paf + aal^3aaluj-36l3aluj^auj^^aujl3l-36all3uj'^ 
-UajPPf + aiaPf + 3aiaLo'^ - Uaico/Spf - 12/3aiLO^, 

Hb = -l2zpa\-2Azpaluj^2zaal + 2zl5d^a\-12alzPuj'^ + a\t 

-Salzf3f3l + Azaajco + 2aitco + Azf3d'^aiuj + 2aiza/3i + 2aizauj'^ 
^SaiLozPPl + 2zPd^uj'^ + Az^lo'^PI ^z + tuj'^ + Az^pf + t^^ + 2zd^PI3l, 



A{t, z) := 768(i^z^aaf/3cj^ - IGti^z^a^a? - 2{)d'^z^a^a\ + Scj^ + GGti^z^at/^ 

-Vimd^z'^aalpuj'^ + Sa^ + 384^i^tza?/3cj^ - md^tzpujai - QAdHzaa\uj^ 

-WdHzaaiu^ + Ilb2d'^z^aalf3uj'^ + 384t/2tza^/3a; + b76dHzajf3uj'^ 

-192dHz(3uj'^al - 192dHzPu^ai - IGd'^z^aaf + 3ald'^ - 2M^t^uj'^al 

-Wdh'^cj^ai - Wd'^t'^alcj - lUd'^z^Paj - AM'^z'^pJ^ + ^M^tzPuo'^al 

-lUd^z^p^aj - lUd^z'^fS^uj^ - SdH^ujai - SdHzaj - SdHzuj^ 

—^M^tzaa\(jj — 9^d^tzaa\uj'^ — 32d^tzaa\uj — Wd'^tzaaiuj'^ + Sizd'^auj'^ 

-mzd^pal + ^izd^aa\ + m^zauj'^ + m^zaal - md^zalp - AdH^uj'^ 

-Sizd'^ai - Sizd'^uo - ASd"^ z^ a'^ afuj - AOd"^ z^ a'^ ajuo'^ - 230Ad'^z^al/3'^uj 

-ASdHzPaj - AMhzPuj'^ + M'^z'^au - Ad^z^a'^cj'^ - Ad^z^o?o!{ - Ad^z'^a'^cj^ + Sd'^z'^aai 

-S6Ad^z^al/3^uj^ - 576d^ z^ai/3^uj^ - 32d^z^aaluj - IGd^z^aaiu^ - AdH^a^ - AdH^u^ 

-AdH^aj - 1152d'^z^al/3^uj - 576d^z^al/3^uj - b76d'^z^a^/3^uj^ + Qd^aiuj - Sd^ z'^ a'^ aiuj 

-96id'^zai/3uj'^ - 192id'^zal/3uj + Wid'^zaaiu - 3SAizd^ /3ajuj - 38Aizd^ pajuj^ 

-\-32izd^aaluj - 576izd^ I3aluj^ + ASizd'^aco^aj + 32izd'^au;^ai - 230 Ad^z^a^ 13^ uj^ 

-3A56d'^z^alP'^(jj'^ - b7M^ z'^a\p'^uj'^ + 3d^uj'^ + I92d^z'^aalp 

-96izd'^ Paiou^ + 12ou^ai + IScefo;^ - GAd'^ z'^ a'^ alou - ASdHz^cu^ 

-96d^z^a^aiuj^ - 6Ad^ z^ afuj^ - 16d^ z^ ajuj"^ + ASdHzafP 

-^192d^z'^aalf3 - 16d^ z'^ a'^ aiu^ + 9MHzalp - IQdHzaal 

-IGdHzaal - 576d^z^ali3^ - 576d^z^al/3^ - lAAd^z^^^uj^ - lAAd^z^a^^^ 

-2SSd^z'^P'^(jjai + 192d'^z^alf3uj + 96d'^z'^alP(jj'^ + Uafuj - l^dHzooai 

-192d^z'^P(jjai + 76Sd'^z'^aalP(jj + 192d^z'^aalPuj'^ + 3Md'^ z'^ aa\Puj - Az'^d^, 
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-192(fz^aal(3uj^ ^ af - 384dhzal(3uj^ + 96dHz(3ujai + 64dHzaaluj^ 
-^IQdHzaaiuj^ - 1152 d'^z^aal (3 uj'^ - 384dHzal(3uj - 576dHzaf(3uj'^ 
^192dHz(3uj'^al + 192dHz(3uj^ai + IGd'^z'^aaf + ajd'^ + 24dH'^uj'^al 
^Wd'^t^uj^ai + Wd^t^aluj + 144d'^z^/3a? + 48(i'^^^/3cj^ - 96d^tzPuj^al 
-^144d^z^P'^al + 144(i^z^/3^cj^ + SdH^ujai + MHzal + MHzuo'^ + 64d^tzaatcj 

^40d^z^a'^aluj'^ + 2304^^2z^a^/3^a; + 48d^tz/3a? + 48dHz(3uj'^ - M^z'^auj 
^4d''z^aW + 4d''z^a^al + 4dS2a2c^^ - Sd^z^aa^ + 864^^^^ a? /^^c^^ 
+576d^z^ai/3^cj^ + 32d^z'^aaluj + IGd^z^aaicj^ + 4(i^t^at + 4d'^t^uj^ 
^4dH'^uj^ + 4^"^^^^^ + 1152d^z^al/3'^uj + 576d^z^a?/3^cj + 576d^z'^af/3^uj'^ 
-^2d^aiuj + 8d^z^a^aicj + 2304d'^ z'^a^P^uj^ + 3456(i^z^a?/3^a;^ 
+576(i^z2at/3^a;'^ + d'^u'^ - 192d^ z'^ aalp + 4cj^ai + 6a?a;^ 
+64^^2z2a2a^a; + 48d'^tz/3a;'^ + 9Qd^ z'^ a\uj'^ + 64^^2z^a2a?cj^ 
+16t/^z^a^a?cj'^ - 4MHza\p - 192d^z^aalP + Wd"^ z'^ aiuj^ 
-9Qd^tzalp + Wd'^tzaal + IddHzaal + 576^i^z^a?/3^ + 576(i'^z^a?/3^ 
+144^/^^2/3^0;'^ + 144(i^z2at/52 + 288(i^z2^2a;ai - 192d'^z^al/3uj 
-96d'^z^alPuj^ + 4a?cj + Wd^tzuoai + 192c^'^z2/3cjai 
-768d2z2aa?/3cj - 192d'^z^aal(3uj'^ - 3Md^ z'^ aa\puj + 4z'^d^. 
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Fig. 1: (color online) One solition solution {E,p, rf) of the H-MB system with cj = 1.5, ai = 0.5, /3i = l,a = 2,/3=— 1. 
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Fig. 3: (color online) One positon solution {E,p, 77) of the H-MB system when = 1.5, ai 





Fig. 5: (color online) Breather solution (E,p,r]) of the H-MB equation when 
oj = 0.5,(i = 0.5,ai = -1, /3i = 1, a = 2, /3 = -1. 
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Fig. 6: (color online) Breather solution {E,p,r}) of the (0,l)-CMKdV-MB equation when 




Fig. 7: (color online) Breather solution 77 for different values of {a, /3) of the H-MB system when 

UJ 0.5,d = 0.5,ai 1. 




Fig. 8: (color online) Rogue wave solution (E,p,ri) of the H-MB system when 
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Fig. 9: (color online) Rogue wave solution {E^p^rf) of the H-MB equation when 
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Fig. 10: (color online) Enlarged graph of rogue wave solution p (in Fig. 9| of the H-MB system when 

oj = 0.5,d = 0.5,ai = = l,a = 2,/3 = 
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Fig. 11: (color online) Rogue wave solution {E^p^rf) of the (0,l)-CMKdV-MB equation when 
oj = 0.5,(i = 0.5,ai = -1, /3i = 1, a = 0, /3 = 1. 
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